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ABSTRACT 

The mass of the central black hole in the giant elliptical galaxy M84 has previously been measured by 
two groups using the same observations of emission-line gas with the Spac e Telescope Imaging Spectrograph 
(STIS) on the Hubble Space Telescope, giving strongly discrepant results: iBower et alj (11998b found M BH = 
(1.5iJ g) x 10 9 M , while iMaciejewski & Binnevl ( 120011) estimated M BH = 4 x 10 s M . In order to resolve this 
discrepancy, we have performed new measurements of the gas kinematics in M84 from the same archival data, 
and carried out comprehensive gas-dynamical modeling for the emission-line disk within ~ 70 pc from the 
nucleus. In comparison with the two previous studies of M84, our analysis includes a more complete treatment 
of the propagation of emission-line profiles through the telescope and STIS optics, as well as inclusion of the 
effects of an intrinsic velocity dispersion in the emission-line disk. We find that an intrinsic velocity dispersion 
is needed in order to match the observed line widths, and we calculate gas-dynamical models both with and 
without a correction for asymmetric drift. Including the effect of asymmetric drift improves the model fit to 
the observed velocity field. Our best-fitting model with asymmetric drift gives M BH = (8.5^'g) x 10 8 M Q (68% 
confidence). This is a factor of ~ 2 smaller than the mass often adopted in studies of the M B e - 0* and M B h —L 
relationships. Our result provides a firmer basis for the inclusion of M84 in the correlations between black hole 
mass and host galaxy properties. 

Subject headings: galaxies: active - galaxies: individual (M84, NGC 4347) - galaxies: kinematics and dynam- 
ics - galaxies: nuclei 



1. INTRODUCTION 

Located in the Virgo cluster, M84 (NGC 4374) is an ellipti- 
cal galaxy with a radio - loud a ctive galactic nucleus (AGN). 
According to iHo et alj dl997h . the AGN is optically clas- 
sified as a Type 2 low-ionization nuclear emission-line re- 
gion (LINER). M84 is k nown to contain a circumnuclear gas 
disk ([Bower et all 1 19971) . providing the opportunity to mea- 
sure the mass of the central black hole if the gas partici- 
pates in circular rotation in a thin disk-like structure. Fol- 
lowing the installation of the Space Telescope Imaging Spec- 
trograph (STIS) on the Hubble Space Telescope (HST), M84 
was the first targ et for a gas-dynam ical measurement of its 
black hole mass ( IBower et alj 1 1998b . Gas-dynamical mea- 
surements had previously been made for several objects using 
data from the Faint Object Spectrograph (FOS) and Faint Ob- 
ject Camera (FOC) dHarms et al.lll994tlFerrarese et alJ | 1996b 
Macchetto et al. 1997; I van der Marel & van den Bo sch 1998; 
Ferrarese & Ford! [19991) . but STIS provided a dramatic im- 
provement in data quality for spatially-resolved emission-line 
s pectroscopy. 

IBower et ail dJ998) showed that the M84 gas kinematics 
could be modeled as a disk in circular rotation. However, near 
the M84 nucleus, the STIS spectra are very complicated with 
severely blended, double-peaked, and asymmetric Ha and 
[N II] line profiles. Bow er et alj {f998) interpreted the com- 
plex nature of the emission lines as arising from two kinemat- 
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ically distinct gas components. One high velocity gas com- 
ponent was thought to be part of the nuclear gas disk rotat- 
ing about the black hole, while the other low velocity compo- 
nent was believed to b e an extension of a filamentary structure 
seen on larger scales. IBower et al.1 d 19981) therefore included 
only the high velocity component in their modeling, and de- 
termined a black hole mass of (1 -5t^\) x 10 9 M . 

An alternative explanation for the complex line pro- 
files seen near the M84 nucleus was suggested by 
Maci eiewski & Binnevl ([2001). They argued that a rapidly 
rotating disk could give rise to the complex line profiles 
and that the complicated spectra are due to a slit width that 
is larger than the core of the telescope point-spread func- 
tion (PSF). The location at which light enters the slit af- 
fects the inferred velocity, creating an instrumental veloc- 
ity gra dient along the dispers i on dir ection. In the case of 
M84, Macieiewski & Binnevl (120011) identified a "caustic" 
that occurs where the instrumental velocity offset balances the 
change in velocity due to the gas rotating about the black hole. 
They used the location of the caustic to estimate a black hole 
mass of 4 x 10 8 M Q . Similarlv. lBarth et al.ld2001bl) attempted 
to model the complex line profiles in M84 and found that the 
double-peaked and asymmetric profiles could be qualitatively 
reproduced with a single, rapidly rotating disk component ob- 
served through the STIS 0"2-wide slit, however they were 
unable to determine a black hole mass because a direct fit of 
the modeled line profiles to the data was not well constrained. 

M84 has a stellar velocity dispersion of 296 ± 14 km s" 1 , 
and the black hole in M84 lies near the high-mass end of the 
correlations between black hole mass and host-galaxy proper- 
ties, such a s those with the bulge stel l ar velocity dispersion 
(M BH - a*; iFerrarese & Merrittl 120001; iGebhardt et all l2000t 
iTremaine et all l2002t iGiiltekin et alj 120091) and bulge lumi- 
nosity (M BH - L: lKormendv & Gebhardtll2001l; IGiiltekin etaT] 
2009). These relationships imply that the growth of black 
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holes and bulges are intimately linked, and the correlations 
have crucial implications for galaxy formation and evolution. 
Properly determining black hole masses at the high-mass end 
of the relationships is particularly important because those 
measurements will have a substantial effect on the slope and 
scatter of the relationships, which in turn affect the inferred 
black hole mass function. Additionally, accurate mass mea- 
surements of black holes at the upper-end of the Mbh - ct* 
relationship are necessary in order to address questions about 
whethe r the Mbh ~ 0+ or M rh - L relationship is more funda- 
mental (lLauer et al.ll2007l) . 

Recent work has suggested that some black hole masses at 
the high-mass end of the Mbh _ c* and Mbh ~L relationships 
may have been underestimated in previous stellar dynamical 
measurements. In the case of M87, which contains one of 
the largest measured black holes, Gebh ardt & Thomai (120091) 
found that including the galaxy's dark matter halo in the 
stellar-dynamical modeling increased the i nferred black hole 
mass b y a factor of - 2 to 6.4 x 10 9 M Q . IShen & Gebhardtl 
(2010) re-examined the stellar kinematics of NGC 4649, an- 
other galaxy harboring a high-mass black hole, and found 
a factor of ~ 2 increase to 4.5 x 10 9 M Q from the original 
stellar-dynamical mass measurement. Unlike M87, the addi- 
tion of a dark matter halo had a minor effect on the black hole 
mass. Instead, the mass increase was attributed to the differ- 
ent orbital sampling - the orbital coverage from the previous 
model did not adequat ely sample the phase space occupied by 
tangential orbits. Also. lvan den Bosch & de ZeeuwU2010l) in- 
vestigated the effects of using triaxial stellar-dynamical mod- 
els for a couple of objects that were originally modeled as 
axisymmetric systems. Using a triaxial geometry in place 
of an axisymmetric shape had a significant effect on one of 
the objects (NGC 3379), and the black hole mass increased 
by a factor of ~ 2 to 4 x 10 8 M Q . While the black hole in 
NGC 3379 does not fall a t the high-mass end of the M B h - o* 
and Mrh —L relationships. Ivan den Bosch & de Zeeuwl (120101) 
note that triaxial systems may be prevalent at the high-end of 
the relationships. 

These issues provide a renewed motivation to pursue gas- 
dynamical measurements in massive, early-type galaxies. 
Since gas-dynamical measurements of black hole masses rely 
on gas in circular orbits at small radii, within the black hole's 
dynamical sphere of influence, they are insensitive to the 
large-scale effects of a dark matter halo or stellar orbital 
anisotropy. Thus, they can serve as an important cross-check 
to the much more complex stellar-d ynamical models. 

The spectra of M84 illustrated bv lBower et all dl998i ) show 
that the black hole's dynamical sphere of influence is ex- 
tremely well resolved in the STIS data; this is demon- 
strated by the nearly Keplerian falloff in rotation velocity 
with distance from the nucleus. This makes M84 a partic- 
ularly valuable object for constraining the upper end of the 
Mbh _ o* relation. The discr epancy between the masses de- 
rived b y iBower et al.l dl998l) and by IMaciejews ki & Binnevl 
d2001l) has remained troubling, and th ere are a var i ety of 
possible explanations. The modeling of IBower et al.l (119981) 
was based on the decomposition of the emission-line pro- 
files into rotating a nd non-rotating compo n ents, but the 
models calculated by IMaciejewski & Binnevl (120011) clearly 
illustrated how a rotating disk observed through a wide 
slit could give rise to the complex li ne profiles seen in 
M84. Also, the models calculated by IBower et all (119981) 
did not include all of the detailed effects relevant to prop- 



agation of emission-line pr ofiles through the STIS optic s 
that were later ex p lored by IMaciejewski & Binnevl d2001l) . 
iBarth et al.l (120011) . iMarconi et all d2003l) . and others; these 
effects could plausibly have a substantial impact on the de- 
rived black hole mass. On the other hand, the mass derived by 
Maciejewski & Binney (2001) was based on the visual identi- 
fication of a caustic feature in the velocity field and not from 
quantitative fits of models to the data, hence the uncertainty 
in their estimate of Mbh is diff icult to gauge - 
Furthermore, nei th er IBower et al.l 0998) nor 
Maciejewski & Binney (2001) included any possible ef- 
fects due to an intrinsic velocity dispersion in the disk 
or asymmetric drift in their modeling. Many nuclear 
gas disks in early-type galaxies exhibit an internal veloc- 
ity dispersion that may be dynamically important (e.g., 
van der Marel & van den Bosch] 119981: IVerd oes Kleij n et al.l 
2002t iNeumayer et al.1 120071) . The physical origin of the 
intrinsic velocity dispersion is not understood, and early 
gas-dynamical models did not include the possible effect of 
the intrinsic velocity dispersion on the black hole mass. In 
contrast, more recent gas-dynamical models account for the 
possibility that the intrinsic velocity dispersion provides dy- 
namical support to the disk dBarth et al.ll200il ICoccato et al.l 
2006; Neu mayer et al.ll2007l) . The extremely broad emission 
lines near the center of the M84 disk suggest that there could 
be a dynamically significant internal velocity dispersion, 
which could have a substantial impact on the measured black 
hole mass. 

In this paper, we aim to resolve the uncertainty in the M84 
black hole mass by revisiting the archival STIS observations 
and carrying out more comprehensive dynamical modeling 
than has previously been attempted for this galaxy. We de- 
scribe the archival HST STIS observations in §|2]and the mea- 
surement of the emission lines in $3] We then present the 
observed emission-line velocity, velocity dispersion, and flux 
in SI In §E we discuss the specifics of our gas-dynamical 
model. We provide details of the stellar mass profile, the 
emission-line flux distribution used in the model, the contri- 
butions from various sources to the line widths, and the asym- 
metric drift correction. The final disk models (both with and 
without an asymmetric drift correction) are presented i n ^6 
and we quantify the possible sources of uncertainty in {36.3 
Finally, in $7] we compare our black hole mass measurement 
to past results. Throughout the paper, we adopt a distance to 
M84 o f 17 .0 Mpc in orde r to be consistent with IBower et al.l 
(fl998h and lGiiltekin etail (120091 

2. OBSERVATIONS AND DATA REDUCTION 

M84 wa s obs erved under program GO-7124 (see 
IBower etal.1 119981) with the STIS 52x0.2 aperture and 
no gaps between the three adjacent slit positions. The slit 
was aligned at a position angle of 104° east of north, approx- 
imately perpendicular to the radio jet and near the major axis 
of the gaseous disk. The CCD was read out in an unbinned 
mode giving a wavelength scale of 0.554 A pixel -1 and a 
spatial scale of 0''0507 pixel -1 . The G750M grating was used 
to provide coverage of 6295-6867 A, which includes the 
Ha, [Nil] AA6548,6583, and [S II] AA6716,6731 emission 
lines. The exposure times ranged from 2245 to 2600 s per slit 
position, with two individual exposures at each position. 

The data were reduced using the standard Space Telescope 
Science Institute (STScI) pipeline. The pipeline includes dark 
and bias subtraction, flat-field corrections, and the combining 
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of subexposures to reject cosmic rays. The data were flux 
and wavelength calibrated, and rectified for geometric distor- 
tions. Before the geometric rectification, we performed an 
additional cleaning step to remove any hot pixels or cosmic 
rays that remained in the flux-calibrated images. 

We also obtained Wide Field Planetary Camera 2 (WFPC2) 
F547M, F658N, and F814W images from the HST archive, 
originally observed under program GO-6094. In each image, 
the nucleus was centered on the PC detector. The images for 
each of the filters were taken as a sequence of two individ- 
ual exposures, and we used the IRAF combine task to av- 
erage the exposures together and reject cosmic rays. How- 
ever, a number of cosmic rays remained even after using the 
combine task, and so we applied an extra cosmic ray clean- 
ing step to the comb ined image using the LA-COSMIC task 
(vanDokkum 2001). The total exposure times were 1200 s, 
2600 s, and 520 s for the F547M, F658N, and F814W images, 
respectively. 

We created a continuum-subtracted Ha+[N II] image by 
subtracting a scaled combination of the F547M and F814W 
images from the F658N image. We experimented with differ- 
ent scaling factors, and searched for a set of factors that would 
produce a background region in the continuum-subtracted im- 
age with a mean flux as close to zero as possible. In FigureQ] 
we show the continuum-subtracted image with the location of 
the STIS slits overlaid, as well as the F547M W FPC2 image. 
As previously described by Bower et al. ( 1997), the emission- 
line image reveals a compact central source surrounded by an 
extended disk-like structure that traces the nuclear dust disk. 

3. MEASUREMENT OF EMISSION LINES 

We extracted spectra from individual rows of the 2D STIS 
image out to about 0"85 from the slit center. The rows were 
extracted as far out as the emission lines were detectable. Far 
from the slit center, where the emission lines were weak, we 
binned together multiple rows to improve the S/N. Before fit- 
ting the emission lines, we removed the continuum from the 
spectrum. For each row, we fit a line to the continuum regions 
between rest-wavelengths 6500-6650 A and 6700-6750 A, 
and then subtracted the continuum fit from the spectrum. Gen- 
erally, a straight line is not a good description of the contin- 
uum near Ha, but with the low S/N in the continuum at most 
positions and the small wavelength range, we could not per- 
form a more accurate subtraction. 

We applied a Levenberg-Marquardt least-squares min- 
imization techn ique using the MPFIT library in IDL 
dMarkw ardt 2009) in order to fit each spectrum with a set 
of Gaussians, foll owing the same technique described by 
Wals hetail d2008b . We used the propagated error spectra 
to weight the data points during the fit. For all but the 
three innermost rows of the central slit, we simultaneously 
fit five Gaussians to the Ha, [N II] AA6548,6583, and [S II] 
AA6716,6731 emission lines. All the lines were required to 
have a common velocity, and the fluxes of the [N II] lines 
were held at a 3 : 1 ratio. The widths of the [N II] lines were 
required to be the same. Also, the [S II] lines were constrained 
to have equal velocity widths. 

Spectra extracted from the central slit position at 0."15 and 
0."10 from the nucleus clearly exhibit double-peaked line pro- 
files. Our basic model described previously provides a good 
fit to the spectra at 0"15 and 0."10, but we also tried fitting two 
Gaussians to each emission line (representing two kinemat- 
ically distinct gas components). While the two-component 
model fits the spectra very well at these locations, the x 2 P er 



degree of freedom (xt) decreased by less than 12% compared 
to the basic model. 

The basic model, however, did not return adequate fits for 
the spectra extracted from the innermost three rows of the 
central slit position, which have a very complex, asymmet- 
ric shape with severely blended Ha and [N II] lines that may 
be double-peaked. We then applied an additional constraint to 
the basic model that required the widths of the Ha and [N II] 
lines to be equal. Even with this model, the fit was poor and it 
was impossible to obtain accurate and unique mean velocities, 
velocity dispersions, or fluxes. The velocity dispersion values 
were very large, ranging from ~ 600 - 900 km s" 1 , and the 
best-fit Gaussian appeared to be systematically shifted blue- 
ward of the peak [N II] A6583 flux in all three rows. We also 
tried fitting a number of more complicated models, such as 
adding a single Gaussian broad component to the Ha + [N II] 
complex or fitting a two-component model. These more com- 
plex, multi-component fits were not well constrained by the 
data, often resulting in Gaussian components with zero flux. 
Moreover, including a broad component in the fit usually re- 
sulted in unreasonably small fluxes and line widths for the 
Ha and [N II] narrow components, although these models did 
seem to provide a better estimate of the mean velocity than 
the previous fitt ing attempts . 

Past work bv lBower et al.l (fl998) has interpreted these com- 
plex line profiles as arising from two separate gas compo- 
nents. They fit two Gaussian components at number of lo- 
cations within 0."3 from the nucleus, and only used one of the 
two components in order to derive a black hole mass. How- 
ever, a s des cribed above, mode ling by Macieiews ki & BinnevI 
(120011) and iBarth et al l (1200 lbl) has shown that the line pro- 
files can be described by a rotating disk model rather than 
requiring two dynamically distinct components. Our detailed 
modeling further supports this interpretation and will be dis- 
cussed in $6] Thus, for the spectra located at 0."15 and 0"10, 
we use the fit results from the basic model, with one Gaussian 
component for each narrow emission line. 

In the three innermost CCD rows of the central slit position, 
the complexity and blending of the line profiles prevents an 
accurate decomposition into Gaussian components, and it is 
not clear whether a broad Ha component is even present. The 
broad wings on the Ha+[N II] complex could plausibly result 
from either a compact broad-line region, or from the rapid 
rotation of the innermost, unresolved portion of the emission- 
line disk, or both. Thus, we consider all measurements from 
these central three rows to be unreliable, and we do not use 
them to constrain our dynamical models. Since the line pro- 
files in these central three rows are dominated by PSF blurring 
and rotational broadening, they add relatively little useful in- 
formation to constrain the disk models, and we will show in 
^6]that the modeling results are insensitive to the values of the 
observed velocity from these central three STIS rows. How- 
ever, for reference in the remaining figures presented in the 
paper, we continue to plot the mean velocity (measured from 
a fit using the basic model along with a single Gaussian broad 
component and a constraint requiring the widths of the narrow 
Ha and [N II] lines to be equal). We do not show the velocity 
dispersion or flux measurements from the central three rows 
in these figures though because we were unable to obtain rea- 
sonable measurements from any of our fitting attempts. 

In Figure [2] we show a few examples of the basic model fit 
to the spectra extracted from the central slit position at ±0"1 
and ±0"15 from the nucleus. We also display the fits to the 
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FIG. 1. — HST WFPC2/PC F547M (left) and continuum-subtracted Hq+[N II] (right) images of the M84 nucleus. STIS observations were made at the three 
adjacent slit positions shown by the rectangular regions overplotted on the continuum-subtracted image. The length of the rectangular regions corresponds to the 
region over which we could measure the emission lines. Both images have been rotated such that the STIS instrumental y-axis points upward. Each box is 3"7 
on a side. 
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FIG. 2. — Spectra extracted from the central slit located within 0." 15 from the M84 nucleus. For the spectra located at ±0" 1 and ±0"15 from the nucleus (left), 
each emission line was fit with a single Gaussian according to the basic model described in iJ5] For the spectra extracted from the three innermost rows (right), 
we show the fit results from the basic model with an additional constraint on the Ha and [N II] line widths, and the same model but with a single Gaussian broad 
component (shown by the blue solid line). In all panels, the single Gaussian narrow components are shown by blue dashed lines and the sum of all the Gaussian 
components are given by the red solid line. 



spectra from the innermost three rows using the basic model 
plus an additional line width constraint, as well as the basic 
model along with a single Gaussian broad component and a 
line width constraint on the Ha and [N II] narrow components. 

In order to estimate the uncertainties on the free parameters, 
we used a Monte Carlo technique. For each row, we generated 
new spectra by adding random Gaussian noise to the original 
spectrum, where the la level of the perturbation was set by 
the residuals between the best-fit model and the original spec- 
trum. The new spectra were then refit, and the 1 a uncertain- 
ties were taken to be the standard deviation of the distribu- 
tion. This method resulted in velocity and velocity dispersion 
errors that were typically 27% and 40% larger, respectively, 
than the formal uncertainties from the best-fit model. 



4. THE OBSERVED VELOCITY FIELD 

From the fits to the [N II] A6583 emission line, we obtained 
the velocity, velocity dispersion, and flux, and we plot these 
values as a function of position along the slit in Figure[3] We 
defined the galaxy center (Y-Offset = 0") to be the row with 
the largest continuum flux, which also coincided with the peak 
of the narrow [N II] A6583 emission-line flux. The grey open 
squares mark the three central rows where the complex line 
profiles precluded any definitive measurements of the [N II] 
mean velocity, velocity dispersion, or flux. These grey points 
mark the velocity of the [N II] emission peak as determined by 
the fits that included a broad Ha component. For reference, 
we also show the best-fit systemic velocity determined from 
our gas-dynamical modeling. 

The multi-slit velocity curves show the gas participates in 



5 



V (km s- 1 ) 



<j([N II]) (km s" 



f([N II]) A6583 
(10~ 16 ergs s _1 cm -2 ) 



1200 


1 , 
- 1 


1 | I 1 1 1 | 1 1 1 I | I 1 1 I |_ 

I — ■ — | ^*at> 

' •** — — — — - 


800 






1200 


- 2 




800 




id 


1200 


1 

" 3 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 

*** 


800 


_ 1 1 1 


, 1 , , , , 1 , , , , 1 , , , , f 



100 



1000 



100 



1000 



100 



r 1 

...-•-•*'••••••• 

: l 1 1 1 1 




-2 

- »* *„ 


1 


=111111111111 


111 = 

■l-H mim_ 

, , , , i: 




■1 -0.5 0.5 1 
Y-Offset (") 



•1 -0.5 0.5 
Y-Offset (") 



1 -0.5 0.5 1 
Y-Offset (") 



FIG. 3 . — The velocity, velocity dispersion, and flux measured from the [N II] A6583 emission line as a function of location along each of the three slit positions 
in M84. The grey open squares mark the velocity measurements from the three central rows where we had difficulty fitting the spectra, while the unreliable 
velocity dispersion and flux measurements from the three innermost rows are left off the figure. The numbers at the top left of each plot correspond to the slit 
locations depicted in Figure [T] The positive and negative Y-Offsets correspond to spectra extracted from the rows above and the rows below the central row, 
respectively, when the slits are oriented such that the STIS instrumental y-axis points upward. Thus, the bottom half of the slit drawn in Figure [T] corresponds 
to negative Y-Offsets and the top half of the slit in Figure [T] corresponds to positive Y-Offsets. Spectral rows that were binned together in order to improve the 
S/N are shown with error bars in the Y-Offset direction. For reference, the dashed line shows the best-fit systemic velocity determined from our gas-dynamical 
modeling. 



regular rotation. There is a steep velocity gradient across the 
inner 0."2 of slit position 2, where the radial velocity drops 
from 1380 km s" 1 at Y-offset = -0."1 to 740 km s" 1 at Y- 
offset = 0." 1 . The radial velocities measured from slit position 
3 show that the gas located to the southwest of the nucleus 
is mostly redshifted relative to the galaxy's systemic veloc- 
ity, while the data from slit position 1 show that the gas to 
the northeast of the nucleus is mostly blueshifted. The [N II] 
A6583 line widths measured from the central slit position are 
large, with velocity dispersions of ~ 400 km s" 1 and ~ 300 
km s -1 at Y-offset = -0."1 and 0."1, respectively. The line 
widths continue to rise toward the nucleus. However, due to 
the complexity of the spectra, which are severely blended to- 
gether due to the large rotational broadening and the PSF blur- 
ring, we were unable to obtain reliable measurements of the 
velocity, velocity dispersion, or flux from the innermost three 
STIS rows. 

5. MODELING THE VELOCITY FIELD 

The velocity field was modeled assuming that the gas is 
in circular rotation in a thin disk-li ke structure. We follow 
a method similar to that outlined by iBarfh'et al. ( 200l]), and 
we refer the reader there for a more detailed discussion of the 
model. At each radius r in the disk, we determine the circular 
velocity relative to the systemic velocity (v sys ), based on the 
enclosed mass M{r), which depends on the black hole mass 
(Mbh), the stellar mass profile, and the stellar mass-to-light 
ratio (T). The stell ar co ntribution to the circular velocities 
will be discussed in §15.11 

We then project the disk velocity field onto the plane of 
the sky for a given value of the disk inclination i in order to 
determine the line-of-sight projection of the rotation veloc- 
ity at each position in the inclined disk. We generated the 



model velocity field on a highly subsampled pixel grid; each 
model pixel was 10 x oversampled relative to the STIS pixel 
size. Previous studie s have used smal l er subsampling fac- 
tors of 5 = 2 and 5 = 4 dBarth et al .11200 It ICoccato et"aLll2006t 
IWold et al.ll2006l) . However, we found that a large subsam- 
pling factor was necessary in order to sufficiently capture the 
considerable changes in velocity that occur over small spatial 
scales near the M84 nucleus, and to more accurately model 
the steep emission-line flux profile. We experimented with a 
range of subsampling factors (s), which will be discussed fur- 
ther in §|6] before settling on a subsampling factor of s = 10. 

We calculate the intrinsic line-of-sight velocity profiles on a 
velocity grid with a bin size that matches the STIS pixel scale 
of 25.2 km s" 1 . The intrinsic line-of-sight velocity profiles are 
assumed to be Gaussian before passing through the telescope 
optics. The Gaussian profiles are centered on the projected 
line-of-sight velocity at each point on the model grid. Fur- 
thermore, the line-of-sight velocity profiles are weight ed by 
the emission-line flux distribution (to be discussed in §I5.21 i. 
The width of th e in trinsic line-of-sight velocity profiles (to 
be discussed in §15. 3b includes contributions from the thermal 
velocity dispersion of the gas, the instrumental line spread 
function (LSF) for STIS, and an intrinsic turbulent velocity 
dispersion. 

The model velocity field is then synthetically "observed" in 
a manner that matches the STIS observations. This synthetic 
observation includes accountin g for the blurring b y the tele- 
scope PSF. We used Tiny Tim (Krist & Hook 2004) to create 
a 10 x oversampled 0"3 x0."3 portion of the full STIS PSF for 
a monochromatic filter passband at 6600 A. Although a 0"3- 
diameter PSF model excludes part of the STIS PSF wings, 
there is a negligible effect on the inferred black hole mass, 
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as will be demonstrated in $6] Each velocity slice of the line 
profile grid is then convolved with the oversampled 0"3 x0."3 
Tiny Tim PSF model. 

After accounting for the telescope PSF, we propagate the 
model velocity field through the STIS slit. The STIS slit is 
allowed to lie at some angle 9 from the projected major axis 
of the gas disk, and the slit can be displaced some distance 
along the slit width (jc ffset) an d along the slit length (y ffset) 
from the black hole. We included the shifts in velocity that 
arise as a result of the finite slit width. The velocity shifts 
occur because a photon's recorded wavelength depends on the 
location along the slit width at which the photo n enters (e.g., 
iBarth et al.ll200"ltlMacieiewski & Binnevll2001l) . 

After the velocity shifts have been made, we rebin the re- 
sulting emission-line profiles to the STIS pixel size. Thus, we 
are left with a model 2D spectral image similar to the STIS 
data. We convolve the model 2D spectral image with the CCD 
charge diffusion kernel. The CCD charge diffusion kernel is 
given by Tiny Tim for model PSFs that are subsampled, and 
it represents the charge that is spread between the immediate 
neighboring (non-subsampled) pixels. 

Finally, spectra are extracted on a row-by-row basis from 
the model STIS 2D image. We fit a single Gaussian to the 
emission line, analogous to the measurements of the emission 
lines from the STIS data. We are thus able to measure the 
model velocity, velocity dispersion, and flux as a function of 
position along the slit. We then determine the best-fit model 
parameters (Mbh, T, i, v sys , 6, Jt ffset, ^offset) that produce a 
model velocity field that most closely matches the observed 
velocity field over a region (r^) that extends from -0."5 to 
0"5 for each slit position. However, we do not include the un- 
certain velocity measurements from the innermost three STIS 
rows in the fit. We measure the black hole mass and the x 2 
from the model fits to the observed velocity only, after sep- 
arately optimizing the model fits to the observed line widths 
and fluxes as described below. The observed [N II] line widths 
and fluxes are not directly dependent on the black hole mass, 
so the black hole mass is determined from the fit to just the 
radial velocity curves. 

5.1. Stellar Mass Profile 

In order to account for the stellar contribution to the grav- 
itational potential, we modeled the intrinsic density distri- 
bution of the central regions of M84 as the sum of spheri- 
cal ly symmetr i c com p onents with Gauss i an profiles, as don e 
by Sar zi et al.l (12001: IBarth et all (l200"lt) . LSarzi et al l d2002h . 
and ICoccato et al.l (120061). Within such a mu lti-Gaussian 
framework (Monnet. Bac on. & Emselleml 119921) . projecting 
the model density profile on the sky to match the observed 
surface brightness is fairly simple given the properties of the 
Gaussian function. This is particularly true when the instru- 
mental PSF is also expressed in terms of Gaussian functions. 
Once the (non-negative) amplitudes of the Gaussian density 
components have been determined, the contribution of each 
Gaussian to the stellar potential can be conveniently evaluated 
in terms of error functions. 

For M 84, we deprojected the V-band surface brightness 
profile of iKormendy et al] (120091) . who combined large scale 
ground-based images with high-resolution HST NICMOS H- 
band data for the central, dusty region. Since the NICMOS 
PSF is fairly complex, we took particular care to describe it in 
terms of Gaussian components by matching to a NIC2 F205 W 
PSF generated with Tiny Tim. A careful representation of the 
PSF is important in the presence of an active nucleus, as in 



the case of M84, because it affects the isolation and removal 
of an AGN component from the stellar mass budget. 

The lumin osity profile of M8 4 contains a compact central 
component. iBower et al.l (120001) show that this central com- 
ponent is due to the AGN itself and not a nuclear star clus- 
ter, so the light from this component should not contribute to 
the stellar mass profile. In order to determine the impact of 
this component on the mass models, we computed two differ- 
ent deprojected stellar mass profiles, one in which this cen- 
tral component was assumed to be AGN light, and another 
in which it was assumed to be starlight with the same stellar 
mass-to-light ratio as the rest of the galaxy. 

In Figure |4j we present the final product of the multi- 
Gaussian deprojection of the surface-brightness profile of 
M84, that is, the stellar contribution to the circular velocity 
in the galaxy. The circular velocity due to a 4 x 10 8 M 
black hole is also illustrated, since this is a plausible lower 
li mit to the black hole mass bas ed on our results and those 
of iMaciejew ski & B inneyl (120011) . The figure illustrates that 
the stellar contribution to the circular velocities is extremely 
small over the radial range probed by our data (~ 0"5), even 
for the lowest reasonable value of Mbh- The stellar contri- 
bution is negligibly small regardless of whether the central 
compact component in the light profile is assumed to be a star 
cluster or an AGN. 

The results shown in Figure [4] indicate that the total mass 
within the inner 0."5 is likely to be dominated by the black 
hole. In other words, the black hole's dynamical sphere of in- 
fluence is extremely well resolved by the STIS observations; 
this is consistent with the nearly Keplerian shape of the ve- 
locity profile along the central slit position. As a result, our 
model fits should be fairly insensitive to the stellar mass-to- 
light ratio T. In fact, as we describe below, our dynamical 
models do not provide any useful constraints on T, and we 
simply assume a fixed value of T for our final model fits. In 
our dynamical models, we used the circular velocities derived 
under the assumption that the central compact component rep- 
resents AGN light. 

5.2. Emission-Line Flux 

As discussed above in |J5] the model line-of-sight veloc- 
ity profiles are weighted by the emission-line flux at each 
point in the disk. Past gas-dynamical models have used 
either an ana lytic form to describe the intrinsic emission- 
line flux (e.g..|yerdoes Kleiin et al.l2002t ICoccato et al.l2006t 
Marconi et al.l l2006t iPastorini et al.l 120071: iHicks & M alkan 

2008) or the continuum-subtracted Ha + [N II] image directly, 
after deconvolution with the telescope PSF (e.g.. IBarth et alJ 
[lOOl tlShapiro et al.[200dlWold et al.l200llDalla BontaTtaLl 

2009) . IBarth et al.l (1200 il) found that folding the deconvolved 
continuum-subtracted image into the model calculations can 
account for small-scale irregularities in the velocity field. 
However, our attempts to use a deconvolved image of M84 
were unsuccessful, and the dynamical model produced poor 
fits to both the observed flux and velocity. It appeared that 
the deconvolved continuum-subtracted image could not ade- 
quately characterize the steepness of the nuclear emission-line 
light on subpixel scales. 

We therefore experimented with analytic forms of various 
complexity for the intrinsic flux distribution. We searched 
for the simplest function that would sufficiently reproduce 
the observed [N II] flux profile while simultaneously pro- 
ducing model velocities closest to the observed velocities. 
We tested parametrizations that were composed of two to 
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FIG. 4. — Stellar contribution to the circular velocity (v c ) obtained from 
the multi-Gaussian deprojection of the M84 V-band surface brightness pro- 
file. The thick black and thin grey radially increasing lines show Vc values 
obtained assuming stellar mass-to-light ratios of T = 4 and 1 , respectively. 
Two different stellar v c values are shown; the solid curves illustrate the case 
where the central compact source in the light profile is assumed to be the 
AGN (in which case it makes no contribution to the stellar mass profile), 
while the dotted curves illustrate the case in which it is assumed to be a com- 
pact stellar nucleus. The black dot-dashed line shows the Keplerian velocities 
produced by a 4 X 10 s Mq black hole, which, when added to the stellar po- 
tential obtained with T = 4, would lead to the red solid and dotted lines that 
trace a nearly parabolic trend. Over the radial range we fit the kinematic data 
(shown by the vertical grey dashed lines, see ^6), v c is almost entirely deter- 
mined by the black hole mass. Whether the light from the central compact 
component is excluded or included in the mass budget does not significantly 
alter the estimate of Mbh ■ 

four components, and the components were either Gaussian 
functions [parametrized as F(r) = i*bexp(-r 2 /2r 2 )], exponen- 
tial functions [parametrized as F (r) = F {) exp(-r/r s )], or con- 
stants. We experimented with models representing intrinsi- 
cally circularly symmetric disks (producing concentric ellip- 
tical isophotes with constant position angle and axis ratio), 
as well as more complicated models where the isophotes of 
the individual components were allowed to have different cen- 
ters, positions angles, and axis ratios. The parameters for each 
possible surface brightness model were determined by com- 
puting disk models including the parametrized surface bright- 
ness model, calculating the resulting model line profiles as 
described above in $5] measuring the emission-line fluxes in 
the resulting model, and then optimizing the fit of the modeled 
fluxes to the observed [N II] flux distribution by minimizing 
X 2 - Ultimately, we found that the best parametrization was the 
sum of four components: three concentric Gaussians describ- 
ing the compact nucleus, plus a more extended exponential 
component with an offset center relative to the positions of 
the Gaussian components. The amplitudes (Fq), scale radii 
(r s ) along the major axis, position angles, axis ratios (b/a), 
and centroid positions for each of these elliptical components 
in the best-fitting model are listed in Table Q] 

5.3. Line Widths 

A number of factors contribute to the width of the line- 
of-sight velocity profiles, such as rotational broadening, in- 
strumental broadening, thermal broadening, and possibly an 



TABLE 1 
Emission-Line Flux Parameters 



Component 


F 


r s 


^cen 


Veen 


PA 


b/a 






(pc) 


(") 


(") 


(°) 




Gaussian 


234.2 


0.5 


0.00 


0.00 


344 


0.9 


Gaussian 


126.6 


2.8 


0.00 


0.00 


33 


0.4 


Gaussian 


6.0 


12.0 


0.00 


0.00 


31 


0.4 


Exponential 


1.0 


50.0 


0.13 


0.02 


69 


0.7 



NOTE. — The amplitude, F, is in arbitrary flux units. The center 
of the ellipse is described by Xceo and _v ui: n, relative to the location of 
the black hole. The position angle, PA, is in units of degrees clock- 
wise, with PA = pointing along the length of the STIS slit. 



intrinsic velocity dispersion in the gaseous disk. Rotational 
broadening occurs because light from different parts of the 
disk falls within the same slit and consequently is blended to- 
gether. Instrumental broadening includes the PSF blurring, 
the velocity shifts that occur as a result of the finite slit width, 
and the effects of charge diffusion between neighboring pix- 
els. All of these instrumental effects, as well as the rotational 
broadening, are explicitly included in our modeling when the 
intrinsic line-of-sight velocity profiles are propagated through 
the telescope and spectrograph. Before the line-of-sight ve- 
locity profiles are propagated through the telescope optics, we 
assign the line profiles a small width resulting fro m the intrin- 
sic instrumental line-spread function estimated by Barth et al. 
(120011) to be ctlsf = 8 km s" 1 and a velocity dispersion of 
<T t h = 10 km s -1 , which is the expected thermal contribution to 
the line width for gas with a temperature of w 10 4 K. Since the 
values of ctlsf and a± are very small compared with the ob- 
served line widths, they have a negligible effect on the model 
calculation. 

Even though the models accounted for the rotational, in- 
strumental, and thermal broadening, our preliminary mod- 
els predicted line widths that were much smaller than the 
observed velocity dispersions in M84. This behavior has 
often been seen in past stud i es of other galaxies (e.g., 
I van der Marel & van den Bosch! Il998t IVerdoes Kleiin et alj 
l200d 120021: IShapiro et all 120061: iDalla Bonta et al.1 120091) . 
However, previous work has also found that rotational and 
instrumental broadening alone can explain the observed 
line widths in some objects (e.g., [M acchetto et al.1 119971: 
Canet tiet all 12005: Atki nson et al.l 12005b Ide Francesco et alJ 
2006, 2008). In order to produce model line widths that match 
the observations, we include a projected intrinsic velocity dis- 
persion of the form 

o-p = cro + CTiexp ■ C 1 ) 

Thus, the width of the line-of-sight velocity profiles is given 
by er m , ctlsf, and a p added in quadrature, combined with 
the width resulting from the propagation of the line profiles 
through the telescope and spectrograph. We determined that 
cto = 77 km s , cti =214 km s" 1 , and ro = 26 pc by fitting 
the line widths predicted from a preliminary model to the ob- 
served line widths. The unreliable velocity dispersion mea- 
surements from the three central STIS rows were excluded 
from the fit. 

The physical nature of the intrinsic velocity dispers ion is 
not understood. Ivan der Marel & van den Bo sch (1998|) pro- 
posed that the intrinsic velocity dispersion is the result of 
local microturbulence, but that the bulk motion of the gas 
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remains in circular motion. Others suggest the local ran- 
dom motions may contribute pressure support to the disk, 
which if ignored will lead to an underestimate of the blac k 
hole mass (e.g., iBarth et aTl 120011: iNeumaver et alJ 120071) . 
IVerdoes Kleijn. van der Marel. & Noel-Storr I (120061) find that 
radio galaxies in particular (including M84) often show evi- 
dence for nuclear velocity dispersions in excess of those ex- 
pected from pure gravitational motions. They attribute the 
excess velocity dispersion to non-gravitational motions in the 
gas. 

While there remains no current clear consensus as to the 
physical mechanism that causes the discrepancy between the 
model line widths and the observed velocity dispersions, we 
examine two scenarios that should cover the possible range of 
masses for the black hole in M84. First, we consider the case 
in which the intrinsic velocity dispersion does not affect the 
circular velocity. Secondly, we analyze the case in which the 
intrinsic velocity dispersion contributes dynamical support to 
the disk, and we apply an asymmetric drift correction. 

5.4. Asymmetric Drift Correction 

If the intrinsic velocity dispersion provides pressure support 
that balances gravity in the gas disk, then the observed mean 
rotation speed (v^) will be smaller than the local circular ve- 
locity (v c ) for a given black hole mass. Models that do not 
account for this effect will lead to black hole mass measure- 
ments that underestimate the true mass. Although it is unclear 
whether the intrinsic velocity dispersion actually does provide 
dynamical support to the disk, we include an asymmetric drift 
correction in our model to estimate an upper limit to the M84 
bl ack hole mass, foll owing the methods previously described 
bv lBarth et all d200lh . 

Briefly, if we assume that the gas motions are isotropic in 
the r and z (cylindrical) coordinates, then asymmetric drift 
correction is given by 

Here, v is the number density of tracer particles or clouds in 
the gaseous disk, and we assume that intrinsic radial velocity 
dispersion a r can be described by an exponential + constant of 
the same form as Equation Q] For a given parametrization of 
ov, we determine <70 and the pr ojected veloc i ty disp ersion fol- 
lowing the method described bv lBarth et al.l (12001b . This pro- 
jected velocity dispersi on is then added in quadrature to <7 t h 
and olsf (discussed in ^5.3b in order to obtain the widths of 
the line-of-sight velocity profiles before propagation through 
the telescope and spectrograph. We found best-fit values of 
(To = 37 km s" 1 , o\ = 255 km s" 1 , and ro = 53 pc by computing 
preliminary disk models with an asymmetric drift correction 
and fitting the model line widths to the observed [N II] line 
widths. A simil ar expression for the asy mmetric drift correc- 
tion is given by IValenzuela et al. I d2007l) . including an addi- 
tional term that accounts for the effect of gas pressure gra- 
dients. We neglect this additional term as it provides only a 
small contribution to the asymmetric drift correction. 

6. MODELING RESULTS 

In order to determine the best-fit model parameters (Mbh, 
T, i, v sys , 9, Xoffset, ^offset)) we minimiz e d y 2 using the down- 
hill simplex algorithm by iPress et all (119921) . Initially, we 



ran preliminary models without an asymmetric drift correc- 
tion allowing all seven model parameters to vary indepen- 
dently. However, it became clear that the stellar mass-to- 
light ratio could not be well constrained by the STIS data. 
The black hole's sphere of influence is so well resolved that 
we are unable to determine a precise value for T. Our 
preliminary models often returned unreasonably low values 
of T = 0. Therefore, w e estimated T ba sed on the B — V 
color for the galaxy using iBell et all (12003b . The B-V color 
from the Third Reference C atalogue of Bright Galaxies (RC3) 
(de Vaucouleur s et al.ll 199 11) is 0.98, which implies T ~ 4 in 
V-band solar units. We then fixed T to this value in our mod- 
els. 

We also calculated preliminary models using different PSF 
sizes, subsampling factors, and fitting regions. During all of 
these model calculations, the parameters Mbh, U 0, v sys , x ffset> 
and yoffset were allowed to vary during the fit, but T was fixed 
at 4 in V-band solar units. Also, the uncertain velocity mea- 
surements from the three central STIS rows were excluded 
from the fit. In Figure [5] we present the results of the models 
without an asymmetric drift correction. 

In order to test the significance of the PSF size on the 
black hole mass measurement, we calculated disk models us- 
ing Tiny Tim PSFs ranging in size from 0."3 to 0"6 in diam- 
eter. All of these model PSFs exclude part of the extended 
STIS PSF wings, with the omitted flux totaling 18%, 13%, 
11%, and 9% of the entire PSF flux for a 0."3, Of 4, Of 5, and 
0."6-diameter PSF, respectively. We found that the PSF size 
has a small impact on the black hole mass measurement, with 
M B h varying between 4.3 x 10 s M Q and 4.5 x 10 8 M for 
the models without an asymmetric drift correction. Because 
the computation time increases rapidly with PSF size, and the 
PSF size has a minimal impact on the black hole mass, we 
used the Tiny Tim 0."3-diameter PSF in our final models. 

When testing the effect of the subsampling factor on the 
black hole mass, we calculated disk models using using sub- 
sampling factors between s = 2 and s = 10 in increments of 
two. We found that the black hole mass did not vary sig- 
nificantly with the subsampling factor. For models with- 
out an asymmetric drift correction, Mbh fluctuated between 
4.1 x 10 s M Q and 4.3 x 10 8 M . Although the black hole 
mass was almost unaffected by the subsampling factor, we 
found that larger subsampling factors resulted in much better 
fits to both the emission-line flux and velocities. Small sub- 
sampling factors were unable to capture the large changes in 
both the velocity and emission-line flux seen near the M84 
nucleus. Thus, in our final models we used s = 10. 

Similarly, we calculated disk models with rg t values be- 
tween 0f'3 and 1". Like with the PSF and subsampling factor, 
the size of the fitting region had a small impact on the black 
hole mass, withM BH between 4.2 x 1O 8 M and 4.5 x 10 8 M o 
for the models without an asymmetric drift correction. Using 
the largest possible fitting region is not necessarily the best 
approach because there can be systematic departures from the 
disk model at large radii, possibly due to disk warping. The 
region in the central slit position near Y-Offset= — 0."5, in par- 
ticular, is not well fit by any of the disk models. We used rg t = 
0."5 in our final models because this region is large enough to 
contain an adequate number of velocity measurements needed 
to constrain the model, but also small enough so that the disk 
model still provides a good description of the observations. 

We also experimented with using different analytic 
parametrizations of the intrinsic emission-line flux distribu- 
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FIG. 5. — Black hole mass as a function of PSF size (with s = 10 and rfj t = 0"5), subsampling factor (with a 0"3-diameter PSF and rfj t = 0"5), and fitting radius 
(with a Of'3-diameter PSF and s = 10), respectively. The quantities Mbh, U 9, Vsys, *offset> an d ^offset were allowed to vary during the fit, but T was fixed at 4 
V-band solar units. All models were calculated without an asymmetric drift correction. 



tion. In Figure [6] we present the results of disk models with- 
out an asymmetric drift correction, in which the emission-line 
flux was modeled with three additional functions of varying 
complexity: 1 Gaussian + 1 exponential (model B), 2 expo- 
nentials (model C), and 2 Gaussians + 1 exponential (model 
D). The isophotes of the components in models B, C, and D 
are elliptical and have different centers, position angles, and 
axis ratios similar to the emission-line flux distribution com- 
posed of 3 Gaussians + 1 exp onen tial adopted for the final 
model (model A) discussed in 35.21 

Models A, B, C, and D reproduced the observed emission- 
line flux within the errors, but the quality of the fit to the ob- 
served velocity varied significantly (xt ranged from 1 1 .4 - 
16.3 for the models without an asymmetric drift correction). 
Although the intrinsic emission-line flux distribution affected 
the quality of the velocity fit, it only had a small influence 
on M B h, with the best-fit mass between 4.2 x 10 8 M© and 
4.7 x 10 8 Mq. The observed emission-line flux could also be 
fit with two additional analytic functions: 3 exponentials and 
3 Gaussians + 1 exponential with a nuclear hole. When op- 
timizing this last model to the observed flux distribution, we 
found a best-fit hole radius of 0.8 pc, which is well within a 
single STIS pixel. The emission-line flux models composed 
of 3 exponentials and 3 Gaussians + 1 exponential with a 
nuclear hole produced a black hole mass of 4.5 x 10 s M Q 
and a xl of 13.7 and 12.3, respectively. These two mod- 
els are left off Figure [6] for clarity, but similarly show that 
the black hole mass doe s not change drastica lly. Our find- 
ings are consistent with Mar coni et al.1 (120061) . who demon- 
strated that the emission-line flux model has little effect on the 
final black hole mass measurement for Centaurus A. Many 
other emission-line flux models were tested, but the models 
returned unacceptable fits to the observed flux, and thus they 
are not presented here. 

The other disk parameters appeared to be well constrained 
by the multi-slit STIS data. Throughout the numerous cal- 
culations of preliminary models without an asymmetric drift 
correction, the parameters i, v sys , 9, JC ffset> and ^offset varied 
at most by 8°, 10 km s" 1 , 6°, 0."03, and 0."03, respectively. 
Similarly, we found deviations of at most 10°, 12 km s" 1 , 8°, 
0"03, and 0."04 for i, v sys , 9, x ffset> and y ffset) respectively, 
for the models with an asymmetric drift correction. The un- 



certainty in the black hole mass due to the se parameters is 
negligible, but will be accounted for in 36. 3 1 when the formal 
model fitting uncertainty is determined. 

6.1. Models Without an Asymmetric Drift Correction 

The final model without an asymmetric drift correction was 
calculated using a 0."3 x 0."3 PSF, a subsampling factor of s = 
10, and a fitting region of r^ t = 0"5. Due to the complexity of 
the spectra extracted from the three central STIS rows and the 
resulting uncertainty in the radial velocity measurements, we 
excluded these three innermost measurements from the final 
fit. The stellar mass-to-light ratio was fixed to a value of 4 in 
V-band solar units. We also inc lude d the projected intrinsic 
velocity dispersion described in §15.31 

In Figure [7] we compare the final model velocity, velocity 
dispersion, and emission-line flux to the observed values as a 
function of location along the slit for each of the three slit po- 
sitions. As can be seen in this figure, the final model is able to 
sufficiently reproduce the shape of the observed velocity field 
and the emission-line flux. It is also apparent that an intrinsic 
velocity dispersion is needed in order to bring the model into 
agreement with the observations. In the left column of Table 
|2] we give the best-fit parameter values for the final model 
without an asymmetric drift correction. 

Even though the final model matches the general shape of 
the velocity curves fairly well, we find that x 2 = 594.7. The 
58 velocity measurements are fit with 6 parameters (Mbh, 
i, v sys , 9, Offset, and ^offset), resulting in \l = 11-4. We ac- 
knowledge that our thin-disk model is unable to reproduce 
all of the velocit y str ucture that is observed. As will be 
shown below in 36.31 removing several velocity measure- 
ments at locations that clearly deviate from pure circular ro- 
tation leaves the black hole mass unaffected while signifi- 
cantly improving xt- Thus, we believe our Mbh measure- 
ment remains credible despite the moderate xl value. Fur- 
thermore, we account for the imperfect disk model in 36.31 
by renormalizing xl before estimating the formal uncertainty 
in Mbh due to the model fitting process. We use this ap- 
proach because it will provide the most conservative esti- 
mate of the uncertainty in the black hole mass. We note 
that several gas-dynamical models of other objects have sim- 
ilar xl values (e.g., ICapetti et alj|2005fc lAtkinson et al.ll2005l : 
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FIG. 6. — Results of disk models with four different parametrizations of the emission-line flux: 3 Gaussians + 1 exponential (model A), 1 Gaussian + 1 
exponential (model B), 2 exponentials (model C), and 2 Gaussians + 1 exponential (model D). We present the model velocity for the central slit position (left) 
and the model emission-line flux from all three slit positions (right), as well as the best-fit Mbh in units of 10 s Mq and the corresponding xi- All models were 
calculated without an asymmetric drift correction using a 0"3-diameter PSF and s = 10. The fitting radius was set to rst = 0"5, and thus the black filled circles 
are the data measured from the STIS observations that were included in the fit while the grey open squares are the data that were excluded. The quantities Mbh > 
i, 6, v S ys, ^offset* and y ff sct were allowed to vary during the fit. 



V (km s" 1 ) 



1200 



800 



1200 



800 



1200 



800 



| I I I I | I I I I | I I I I | 



I I I I I 



I I I I I 



1 1 1 1 1 



1 1 1 1 1 



1000 



100 



1 -0.5 0.5 
Y Offset (arcsee) 



ct([N II]) (km s" 1 ) 



1000 = 



100 - 



000 



100 = 



f([N II]) X6583 
(10~ 1B ergs s _1 cm~ s ) 




■1 -0.5 0.5 
Y Offset (arcsec 



1 -0.5 0.5 
Y Offset (arcsec) 



FIG. 7. — Results of the final model without an asymmetric drift correction using a 0" 3-diameter PSF, s = 10, and ra t = 0"5. The model velocity, velocity 
dispersion, and emission-line flux are shown by the red solid line. The red dotted line shows the model velocity dispersion if an intrinsic velocity dispersion is 
not included in the model. We applied a single scaling factor to the model emission-line flux such that the median values of the model and the data would match. 
The black filled circles are the data from the STIS observations that were included in the fit, while the grey open squares are the data that were excluded from the 
fit. For reference, the horizontal dashed line shows the best-fit systemic velocity determined from the final model. 
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TABLE 2 
Final Disk Model Parameters 



w/o Asymmetric w/ Asymmetric 
Drift Drift 





(4.3^-8) XlO 8 


(8.5+°-9)xl0 8 


T (V-band solar) 


4 


4 


in 




72+' 


v sys (km s" 1 ) 


1060^ 


1060+^ 


en 


27±3 


28+3 


-^offset (") 


02 +0 02 


u - ul -o.oi 


^offset (") 


_ 4+o.02 

-0.01 


_ 05+6.03 



NOTE. — Uncertainties given for the black hole mass are 68% 
confidence limits. The stellar mass-to-light ratio, T, was frozen 
at 4 (V-band solar units) for both models. A relative angle of 
6 = 27° and = 28° corresponds to a disk major axis position 
angle of 77° and 76° east of north. 

Ide Francesco et alj|200llballa Bonta et1ul2009l) . 

6.2. Models Including an Asymmetric Drift Correction 

We explored the possibility that the intrinsic velocity dis- 
persion contributes dynamical support to the gas disk by in- 
cluding an asymmetric drift correction in the model. Similar 
to the final model without an asymmetric drift correction, we 
calculated the disk model using a 0"3-diameter PSF, s = 10, 
and r+it = 0."5. We fixed T to 4 (V-band solar units), and 
excluded the three uncertain central velocity measurements 
from the fit. We also included an i ntrin sic radial velocity dis- 
persion profile a r as described in $5.41 

The asymmetric drift correction depends on the radial gra- 
dient in the number density of clouds in the disk, v(r). We 
took v to be given by the intrinsic emis sion-line flux dis- 
tribution, as has done in the past (e.g., iBarth et aLI 120011 : 
ICoccato et alll2006t iNeumayer et al.ll2007l) . In order to com- 
pute an azimuthally-averaged value for the radial derivative of 
v in the asymmetric drift correction, we fit a revised model to 
the emission-line surface brightness profile, modeling it as an 
intrinsically circularly symmetric, inclined disk. The model 
for v(r) contained three Gaussians and one exponential com- 
ponent (similar to that described in Table[T|i but with the modi- 
fication that all four components were forced to have the same 
centroid, ellipticity, and position angle. Using this model for 
v(r), the asymmetric drift correction v 2 c - was then com- 
puted as a function of radius in each model calculation. The 
disk velocity field is then given by v^. 

When an asymmetric drift correction was included in the 
disk model calculations, the black hole mass increased to 
8.5 x 10 8 Mq. This model is a better fit to the observed ve- 
locities, with x 2 = 524.8 and xt = 10.1, where 58 velocity 
measurements were fit with 6 free parameters. In Figure [8] 
we present the velocity, velocity dispersion, and flux from the 
final model with an asymmetric drift correction along with the 
observed velocity field. We also show the observed 2D STIS 
spectrum and the synthetic spectrum from the final model 
with an asymmetric drift correction for all three slit positions 
in Figure [9] This figure demonstrates that the rotating disk 
model is able to qualitatively match the complicated nature 
of the observed spectra. In the right column of Table ff] we 
provide the best-fit values of the model parameters. 

6.3. Error Budget 

We incorporated a number of sources of uncertainty in order 
to determine the final range of possible black hole masses. In 



addition to the formal model fitting uncertainties, we included 
the uncertainties associated with the stellar mass-to-light ra- 
tio and density profile, the PSF size, the subsampling factor, 
the fitting region, the analytical form of the emission-line flux 
profile, and the velocity measurements from the three central 
STIS rows. We additionally explored the variation in the black 
hole mass when several velocity measurements that appear to 
deviate from circular rotation are excluded from the fit. For 
the disk model with an asymmetric drift correction, we also 
included the contribution to the uncertainty resulting from the 
parametrization of the v. These sources of uncertainty are 
summarized below. 

Model Fitting Uncertainty: We began by estimating the 
model fitting uncertainty. We explored the black hole mass 
parameter space between 3.2 x 10 8 M and 6.2 x 10 8 M Q 
for the disk model without an asymmetric drift correction, 
and from (7.0-9.8) x 10 8 M Q for the model with the cor- 
rection. We held Mbh fixed while the remaining parame- 
ters (i, v sys , 6, ^offset, and y ff set ) were allowed to vary. Be- 
fore calculating these disk models, we rescaled the observed 
velocity uncertainties, which has been the common practice 
in previous work (e.g., ISarzi et all 1200 U IBarth et alj 120011: 
Ide Francesco"et~aT] 120061 2008). We added 29 km s" 1 in 
quadrature to each of the observed velocity uncertainties in 
order to obtain xt ~ 1 f° r the final model without an asym- 
metric drift correction. Likewise, we acquired xt ~ 1 f° r tne 
final model with an asymmetric drift correction by adding 28 
km s" 1 in quadrature to the observed velocity uncertainties. 
By rescaling in this manner, the final uncertainty on the black 
hole mass will increase. This approach provides a conser- 
vative way in which to account for the detailed, small-scale 
velocity structure that our disk models fail to reproduce. 

Figure [10] displays the results of the disk models. The 
range of black hole masses which caused x 2 to increase by 
1.0 from the minimum value provides the 68.3% confidence 
limits (lcr uncertainties) on Mbh- For the disk models with- 
out an asymmetric drift correction, the la uncertainty cor- 
responds to (4.0-4.8) x 10 s Mq, while the range of black 
hole masses is (8.0-8.8) x 10 8 M for the disk models with 
the correction. The 99.7% confidence limits (3er uncertain- 
ties) on Mbh were found by searching for the range of black 
hole masses that caused the minimum x 2 to increase by 9.0. 
The 3cr uncertainties corresponded to (3.5-5.7) x 10 8 M for 
the disk models without an asymmetric drift correction and to 
(7.4-9.6) x 10 8 M Q for the models with the correction. These 
estimates also account for the small uncertainties associated 
with the parameters i, v sys , 9, *offset> and ^offset- 

We found that rescaling the uncertainties on the observed 
velocities resulted in a slightly different best-fit black hole 
mass than the mass found from our fi nal m odel with an asym- 
metric drift correction presented in {36.21 However, the black 
hole mass found from our final model with an asymmetric 
drift correction still falls within the lcr model fitting uncer- 
tainty. Similar behavior has been see n in the previous anal - 
ysis of NGC 1300 and NGC 2748 by lAtkinson etai] (120051) . 
We choose to apply the fractional model fitting error found 
abov e to the best-fit Mbh from the final model presented in 
36.21 Below, we describe the additional contributions to the 
error budget. 

Stellar Mass-to-Light Ratio and Density Profile: Fixing 
T = 4 does not have a large effect on the black hole mass, fur- 
ther emphasizing that the STIS data is fairly insensitive to the 
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FIG. 8. — Results of the final model with an asymmetric drift correction; see Figure[7jfor description. 




FIG. 9. — Continuum-subtracted observed 2D STIS spectrum for the Ha + [N II] region (left) and the synthetic 2D spectrum of the [N II] A6583 A emission 
line from the final disk model with an asymmetric drift correction (right) for each of the three slit positions. The spatial axis is vertical and the wavelength axis 
is horizontal, with increasing wavelengths to the right. Each box is 93 A in the dispersion direction and 2"2 in the spatial direction. The numbers correspond to 
the slit locations depicted in FigurefT] 



stellar contribution. We find that by allowing T to vary during 
the fit, the best-fit stellar mass-to-light ratio is T = and the 
black hole mass is M B h = 4.8 x 10 8 M© for the model without 
an asymmetric drift correction. Thus, increasing the mass-to- 
light ratio from the best-fit value of T = to T = 4 results 
in a 12% decrease in Mbh- For the model with an asymmet- 
ric drift correction, allowing T to vary during the fit caused 
Mbh to change by only 1% from the best-fit mass, while T 
decreased to 2.6 V-band solar units. 

Also, as discussed in 35.11 and shown in Figure |4] the un- 
certainty related to the exclusion of the central compact com- 
ponent from the mass budget is unlikely to affect our Mbh 
measurement. In fact, even if all the nuclear light is assigned 
a mass meaning, the black hole mass remains the same at 
8.5 x 10 8 M and drops to 4.1 x 10 s M (a 5% change from 



the best-fit mass), with and without asymmetric drift correc- 
tion, respectively. 

PSF size: When PSFs of different sizes (Of 3 - 0"6 in di- 
ameter) were used in the calculation of the model without an 
asymmetric drift correction, the black hole mass varied be- 
tween (4.3-4.5) x 10 8 M Q with a root-mean-square (rms) 
scatter of 8.7 x 10 6 M Q . This corresponds to 2% of the best- 
fit black hole mass. Similarly, for the model with an asym- 
metric drift correction, the black hole mass varied between 
8.5 x 10 8 M Q and 8.6 x 10 8 M Q with an rms scatter equal to 
4.3 x 10 6 Mq, or 1% of the best-fit black hole mass. 

Subsampling Factor. For disk models without an asym- 
metric drift correction, using subsampling factors ranging 
from 5 = 2 to 5=10 resulted in black hole masses between 
M B h = 4.1 x 10 8 M© and 4.3 x 10 8 M with a rms scatter of 
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FIG. 10. — Results of models without an asymmetric drift correction (left) and including the correction (right) using a 0"3-diameter PSF, s = 10, and ra, = 0"5. 
For each model, A/bh was held fixed while i, 6, v S y S , ^offset. an d Vofftet were allowed to vary such that x 2 was minimized. The dashed and dotted horizontal lines 
denote where \ 2 has increased by 1.0, and 9.0, from the minimum value, corresponding to la, and 3c, uncertainties on Mbh, respectively. 



only 7.5 x 10 6 M , or 2% of the best-fit black hole mass. The 
black hole mass ranged from (8.2-8.6) x 10 8 M for disk 
models with an asymmetric drift correction. The rms scatter 
was 1.4 x 1O 7 M , which is 2% of the best-fit black hole mass. 

Fitting Region: For the disk models without an asym- 
metric drift correction, the black hole mass varied between 
M BH = 4.2 x 10 8 M Q and 4.5 x 10 s M for fitting regions 
between 0."3 and 1". The rms scatter in Mbh was small at 
1.3 x 10 7 M , or 3% of the best-fit black hole mass. The 
best-fit black hole mass for models with an asymmetric drift 
correction varied between (8.2-8.7) x 10 8 M , with an rms 
scatter that amounted to 1.6 x 10 7 M , corresponding to 2% 
of the best-fit mass. 

Emission-Line Flux: The dynamically cold, thin-disk mod- 
els with different parametrizations of the emission-line flux 
(previously described in © returned best-fit masses between 
4.2 x 10 8 M and 4.7 x 10 8 M . The rms scatter was 
1.6 x 10 7 M , or 4% of the best-fit black hole mass. When 
the different emission-line flux models were used to weight 
the line-of-sight velocity profiles in the model with an asym- 
metric drift correction, we found that the black hole mass fluc- 
tuated between 8.1 x 10 8 M and 8.9 x 10 8 M with an rms 
scatter of 2.7 x 10 7 M , or 3% of the best-fit M BH - 

Velocity Measurements Within 0'.'05 of the Nucleus: We cal- 
culated a disk model without an asymmetric drift correction 
in which the three innermost velocity measurements were in- 
cluded in the fit, as a comparison to the final model that ex- 
cludes these uncertain measurements. The black hole mass 
decreased by just 2% to M B h = 4.2 x 10 8 M . The inclusion 
of the velocity measurements from the three innermost STIS 
rows also had an insignificant impact on Mbh for the model 
with asymmetric drift correction; Mbh decreased by 1 % from 
the best-fit mass. Because the sphere of influence is well re- 
solved by the STIS data, excluding the three central velocity 
measurements during the fit has a small effect on the black 
hole mass. 

Excluding Several Velocity Measurements that Deviate from 
Circular Rotation: The final models with and without an 
asymmetric drift correction match the overall shape of the ob- 
served velocity curves well, but are unable to reproduce all 
of the velocity features, resulting in moderate xt values of 



10.1 and 1 1.4, respectively. In order to measure the impact of 
several velocity measurements that show localized departures 
from pure circular rotation on the black hole mass and xt> we 
fit models with and without an asymmetric drift correction to 
various data sets. We constructed data sets in which we re- 
moved anywhere from a single velocity measurement to ten 
velocity data points that a ppea red t o be discrepant from the 
final models presented in 36.11 and 36.21 such as the points at 
Y-Offset -0.45" to -0.35" from slit position 1, the point at Y- 
Offset -0.55" from slit position 2, and the points at Y-Offset 
0.30" to 0.45" from slit position 3. For the model without 
an asymmetric drift correction, we found that the black hole 
mass varied between 4.1 x 1O 8 M and 4.3 x 1O 8 M , with an 
rms deviation of just 2% of the best-fit mass, and the lowest 
x\ was 6.1. Similarly, for the model with an asymmetric drift 
correction, the rms deviation in the black hole mass was 1% 
of the best-fit mass, varying between (8.3-8.5) x 10 8 M , 
with the lowest xt being 5.3. Thus, the quality of the fit to 
the observed velocity curves improves significantly, while the 
black hole mass remains roughly constant. This lends further 
support that the best-fit masses from our final models with and 
without an asymmetric drift correction are trustworthy despite 
the moderate xt values. 

Number Density of Clouds: For models with an asymmetric 
drift correction, we also tested a variety analytical functions 
for v. These distributions included both shallow and steep 
profiles, and ranged in complexity with anywhere from two 
to six components. The functions were composed of intrinsic 
circularly-symmetric Gaussians and exponentials. When cal- 
culating the models, we continued to weight the line-of-sight 
velocity profiles by the 3 Gaussians + 1 expo nential emission- 
line flux model (model A) discussed in 35.21 The effect on the 
black hole mass was small, and we measured masses between 
(7.1-8.6) x 10 8 Mq with an rms scatter of 5.6 x 10 7 M , 
which is 7% of the best-fit black hole mass. 

All of the above sources of uncertainty were added in 
quadrature to the model-fitting uncertainty. The final range of 
black hole masses for the cold thin-disk model is (3 .6 — 5 . 1 ) x 
10 8 M (la uncertainties) and (3.3-5.8) x 10 8 M (3cr un- 
certainties) with a best-fit black hole mass of 4.3 x 10 8 M . 
For the disk model with an asymmetric drift correction, the fi- 
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nal range of masses is (7.7 — 9.4) x 1O 8 M (la uncertainties) 
and (7.3- 10.0) x 10 s M (3a uncertainties) with a best-fit 
black hole mass of 8 .5 x 1 s M . 

7. DISCUSSION 

By modeling the M84 emission-line gas kinematics as a dy- 
namically cold, thin disk in circular rotation, we find a best- 
fit black hole mass of (4.3±o 8 ) x 10 8 M Q . Incorporating an 
asymmetric drift correction in our disk model results in a best- 
fit black hole mass of (8.5^;|) x 10 8 M . We favor the disk 
model with an asymmetric drift correction because it is phys- 
ically plausible that the intrinsic turbulence affects the disk's 
dynamics and also because the model provides a better fit to 
the observed radial velocities. 

There are a few caveats associated with the asymmetric drift 
correction worth mentioning. Our treatment of the asymmet- 
ric drift correction is only an approximation, and the black 
hole mass measurement becomes increasingly less certain as 
a r /v c approaches ~ 1. In the case of M84, we find that a r /v c 
can reach moderate values of 0.57 (see Figure fTTb. Moreover, 
the increase in the black hole mass when an asymmetric drift 
correction is applied is extremely large (a 98% change), and 
can be attributed to the very steep radial gradients in both v 
and a r . While a r is known fairly well through a fit to the ob- 
served line widths, the radial distribution of the number den- 
sity of clouds in the disk is not known. Following past work, 
we have used the emission-line flux profile as a proxy for the 
density profile v(r). While this is at best a rough approxi- 
mation, we found that different parametrizations of v did not 
result in substantial changes to the inferred black hole mass. 
However, if instead a very simplistic asymmetric drift cor- 
rection is applied, where vj? - vj ot = of, the black hole mass 
decreases by ~ 30% to M BH = 6.0 x 10 8 M . Finally, the 
asymmetric drift correction is strictly applicable in the limit of 
collisionless particles, as in stellar dynamics. Here, we are ap- 
plying the correction to gas clouds, which are not collisionless 
particles, and therefore the analogy is not perfect, although i t 
is still a useful approximation (e.g.. lValenzuela et al. 1 2007). 
Therefore, calculating a model which includes an asymmetric 
drift correction is an informative exercise, but only provides 
an approximate indication of the dynamical influence of the 
intrinsic velocity dispersion. 

Additionally, dust obscuration is an issue in M84 and af- 
fects both our gas dynamical models. As can be seen in the 
WFPC2 V-band image in Figure Q] there are two prominent 
dust lanes within the central ~ 4" that are oriented roughly 
east to west. Based on a V -I color map, lBoweretalJ dl997) 
measure a dust mass of 7 x 10 4 M , and note that the dust ex- 
tinction is the greatest along the northern portion of the central 
dust lane and at a patch located on the nucleus. The effects of 
dust on the determin ation of the stell a r mas s profile should 
be minimized, as the Korme ndy et al.l (|2009) surface bright- 
ness profile uses NICMOS //-band data for the nuclear region. 
Also, as we have shown above, the model fits are insensitive 
to the stellar mass contribution because the black hole dom- 
inates over the region we are modeling. Dust extinction will 
affect the observed gas kinematics, however including these 
effects in the calculation requires radiative transfer models of 
the disk's gas dynamics. Such modeling has yet to be applied 
to gas kinematical measurements and is beyond the extent of 
this paper. 

The best-fit black hole mass of (8.5!g|) x 10 8 M w is a fac- 
tor of ~ 2 smaller than the mass measured by iBower et al.l 
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( 1998) of (1 .5^o g) x 10 9 M Q , although the two measurements 
are formally consistent within their quoted uncertainties. If 
we consider our best-fitting model that does not include asym- 
metric drift (perhaps a more fair comparison since Bower et 
al. did not include asymmetric drift), our result is a factor of 
~ 3 smaller than theirs. While there are many factors that 
contribute to this difference, one reason is likely to be the dif- 
ferent assumptions we have made about the structure of the 
velocity field. Bower et al. decomposed the emission lines 
into two distinct kinematic components, which they character- 
ized as one rotating disk component and one unrelated, low- 
velocity component, while we have modeled the kinematics 
as arising from a single rotating disk. The observed and mod- 
eled line profiles of the central slit position, shown in Figure|9l 
help to illustrate the situation. Close to the nucleus, the emis- 
sion lines flare out into a broad "fan" of emission with large 
line widths and complex, strongly non-Gaussian profiles, due 
to a combination of instrumental and rotational broadening, 
and the incr ease in the intrinsic velocity d ispersion close to 
the nucleus. Maciejewski & Binneyj (120011) describe in detail 
how this fan of emission, which can mimic the appearance of 
two separate kinematic components, can be generated through 
the combination of rotational broadening and the effects of 
observing through a slit that is wider than the PSF core. 

Moreover, our model calculation includes a more detailed 
treatment of t he telescope and sp ectrograph optics than was 
carried out by lBower et al.l (119981) . and includes the full prop- 
agation of emission-line profiles through the spectrograph, 
rather than a simpler propagation of mean velocities. While it 
is difficult to pinpoint the exact cause of the difference in the 
resulting black hole mass, these are likely to be the primary 
reasons. Our results show that we can obtain a robust fit of 
a disk model to the data even in the absence of any measure- 
ments from the central three rows, where the line profiles are 
most severely blended and asymmetric and the velocity field 
is most strongly affected by PSF blurring. 

Of the remaining disk model parameters, our best-fit sys- 
temic velocity and inclination angl e of v svs = 1060 km s " 1 and 
i = 72° differ from the results of Bo wer et all ([1998), who 
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found v sys =1125 km s" 1 and i = 80°. The systemic veloc- 
ity determined through our modeling is in better agreement 
with the recession velocities given by the NASA Extragalac- 
tic Database (NED). The velocities quoted by NED were mea- 
sured from optical lines, and ranged in value from 954 - 1 1 19 
km s -1 , with an average of 1025 km s" 1 . iBower et al.1 (fT998l) 
did not fit for the relative angle between the STIS slit and the 
major axis of the gas disk. Rather, they took the disk position 
angle to be that of the major axis of the large-scal e emission- 
line st ructure (83° east of north) measured by Bau m et ail 
(fl998h . In our disk models, we allowed 9 to be a free param- 
eter and found a similar position angle for the major axis of 
the gaseous disk. Our best-fit value is 8 = 28°, corresponding 
a disk major axis p osition angle of 76° east of north. Also, 
IBower et al.l (1 1998b adopted a stellar mass-to-light ratio of 
T = 5 (V-band solar units), comparable to the one fixed in 
our model of T = 4 (V-band solar units). 

The comparison o f ou r results with those of 
Macieiewski & Binnev (2001) is less direct, since they 
did not carry out model fitting to the observed velocity field. 
Instead, they estimated the black hole mass based on the 
visual location of a caustic feature in the spectrum of the 
central slit position. This caustic results from the interplay 
between the disk rotation and the gradient in instrumental 
wavelength shifts across the slit. These two line-broadening 
effects can be oppositely directed and can effectively cancel 
each other out, resulting in a sharp reduction in the observed 
line width at a specific location in the disk. We do not clearly 
see evidence for this caustic feature in the 2D spectrum 
or in the measured emission-line widths, so we cannot 
reproduce their estimate of Mbh- While such caustics should 
in principle occur in STIS spectra of emission-line disks, it 
is possible that the intrinsic line width in the M84 disk is 
so large that it washes out any observable signature of the 
caustic. We note, however, that our best-fitting model without 
asymmetric drift gives a resul t that is consistent with that of 
iMaciejewski & Binneyl (12001b . 

With respect to the Mbh - 0* and Mbh - L relationships, 
our best-fit black hole mass for M84 of (8.5 ^) x 10 8 M w 
is a fa ctor of ~ 2 smaller than the mass from lBower et al.l 
( 1998) that has been most often adopted in the relationships. 
This new mass measurement lies closer to the mass expected 
from the Mbh — <7* a nd Mbh _ L relationships measured by 
Gulte kin et al.l (12009b . Recent stellar dynamical studies of 
other galaxies, as well as our analysis of M84, show that 
some previous black hole mass measurements should be re- 
evaluated, and determination of the true shape and scatter of 
the Mbh _ and Mbh _ L relationships is an evolving pro- 
cess that will continue to change as observations and model- 
ing techniques continue to improve. Recent calibrations of 
the relations still include some of the early gas-dynamical 
measurements that were done with the FOS on HST, and 
these should be particularly important targets to revisit with 
future observations. Equally important are cross-checks be- 
tween gas-dynamical and stellar-dynamical techniques within 
the same object, as each mass measurement method suffers 
from independent systematic uncertainties. In particular, gas- 
dynamical models would benefit tremendously from a better 



understanding of the physical origin and the dynamical influ- 
ence of the intrinsic velocity dispersion. 

8. CONCLUSIONS 

With the goal of resolving the uncertainty in the M84 black 
hole mass, we have re-analyzed multi-slit archival HST STIS 
observations of the nuclear region of M84. We mapped out 
the 2-dimensional [N II] A6583 emission-line velocity, veloc- 
ity dispersion, and flux. We then modeled the velocity field 
as a cold, thin disk in circular rotation. The line widths pre- 
dicted by this model are smaller than the observed line widths, 
and we found that an intrinsic velocity dispersion is needed 
in order to match the observations. The additional velocity 
dispersion may have a dynamical origin and provide pressure 
support to the gaseous disk, thus we calculated a second disk 
model with an asymmetric drift correction. We found that the 
disk model with an asymmetric drift correction is a better fit 
to the data than the cold thin-disk model. This model gives a 
black hole mass of (8.5^|) x 10 8 M Q (Ict uncertainties). 

We have employed a more rigorous and comprehensive gas- 
dynamical model to measure the black hole mass in M84 
than the previous two studies. In addition to calculating an 
asymmetric drift correction, we included a more sophisticated 
treatment of the effects of the telescope optics and interpreted 
the complex nuclear spectra as arising from a single gas com- 
ponent. We also performed an error analysis that encompasses 
a number of sources of uncertainty in order to determine the 
range of possible black hole masses. 

Our new mass m easurement f o r M8 4 is a factor of ~ 
2 smaller than the IBower et all dl998b value of M B h = 
(1-5+Q g) x 10 9 M Q and a factor of ~ 2 larger than the 
IMaciejewski & Binnevl(l2001b e stimate ofM B H = 4x 10 8 M Q . 
The black hole mass given by IBower et al.l (11998b is often 
used in the Mbh-c* and M^-L relationships. Therefore, 
while recent stellar-dynamical work has found that several 
masses at the high-end of the black hole mass-host galaxy 
relationships have been underestimated by a factor of ~ 2, we 
find that the gas-dynamical black hole mass measurement for 
M84 has been overestimated by a factor of ~ 2. With this ad- 
justment to the black hole mass, M84 lies closer to the mass 
expected from both the Mbh-c* and M^-L relationships. 
Future work should thus aim to understand the systematics 
associated with each of the main mass measurement meth- 
ods, and to re-examine those objects for which uncertain mass 
measurements remain. 
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